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>■ I We discuss the equation of state for cold, dense quark matter in perturbation theory, and 
^ ■ how it might match onto that of hadronic matter. Certain choices of the renormalization 
O ■ scale correspond to a strongly first order chiral transition, and may generate a new class of 
O ! small and very dense quark stars. The results for the mass-radius relation are compatible 
^ I with the recent determination of the mass and the radius of an isolated neutron star by 
O ' Pons et al. 
^ ■ 

^ ! The investigation of the structure of the phase diagram for strongly interacting matter 

pH I is a fascinating subject. Moreover, the study of regions of the diagram associated with 
' extreme conditions in temperature and density can reveal new phenomena, and provide 
a better understanding of Quantum Chromodynamics (QCD). Nevertheless, in spite of 
the increasing interest this field has attracted during the last decade, the mapping of 
the QCD phase diagram is still in its infancy, and most of its regions remain mystifying 
and challenging. An exception is given by the region of nonzero temperature and very 
small densities. There, guidance is given by lattice simulations and experiments in 
heavy ion collisions . On the other hand, an increasing amount of precise astrophysical 
data provides a new resource to probe QCD at large density. The natural candidates to 
study are neutron stars, whose interior might be dense enough to allow for the presence of 
chirally symmetric quark matter Concretely, certain combinations of observables, 

such as the mass-radius relation, may be useful to constrain the equation of state for 
strongly interacting matter. 

Neutron star models predict a maximum mass in the 1.4 - 2.OM0, with a radius 

~ 10 - 15 km. The values for the maximum mass are consistent with the pulsar data 0. 
The conventional theoretical approach used to describe a hybrid star, a neutron star with 
quark matter in the core, or a strange star is the bag model, with at best a correction 
~ Us from perturbative QCD [||. If the bag constant is fit from hadronic phenomenology, 
then the gross features of those stars are very similar to those expected for neutron stars, 
especially in the region of mass compatible with the pulsar observations. Through the 
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years, the study of compact hybrid stars has produced a myriad of different models, with 
all sorts of possible layers and condensates, but with numbers very similar to the ones 
presented above 0. Nevertheless, recent astrophysical results for the mass-radius diagram 
seem to point in a direction unforeseen by traditional models ||^. 

In [|lOl we considered cold, dense quark matter in perturbation theory to ~ @,|], 
analyzing its implications for the physics of quark stars. We used perturbation theory not 
because it is a good approximation, but simply to model the possible nonideality in the 
equation of state of cold, dense QCD. Here we extend this discussion, and consider as well 
how to match this "perturbative" equation of state onto that of hadrons. We emphasize 
especially the comparison between the pressure of the hadronic and quark phases. While 
our results are admittedly similar to other theoretical approaches |]1T|, we try to place 
these results in a more general context, in terms of the QCD phase diagram. Lastly, most 
intriguing is the recent determination of the mass and the radius of an isolated neutron 
star, which suggests an "implausible" equation of state 0. 

Assume that the chiral phase transition occurs at a chemical potential fi^. Our pertur- 
bative equation of state applies in the chirally symmetric phase, when the quark chemical 
potential fi > fi^. In this phase, the effects of a strange quark mass, ^ 100 MeV [^], 
are small relative to the quark chemical potentials, fi > 300 MeV. Thus we take three 
flavors of massless quarks with equal chemical potentials 0|l3l- How the equation of state 
for quarks matches onto that of hadrons is discussed later. 

The thermodynamic potential of a plasma of massless quarks and gluons was calculated 
perturbatively to ~ by Freedman and McLerran |^ and by Baluni [^], using the 
momentum-space subtraction (MOM) scheme. The MOM coupling is related to that in 
the modified minimal subtraction scheme, MS, as [IT0|,O 
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as = g^/{4:7i), with g the QCD couphng constant, and A = 151/48 - (5/18)A^/, with Nf 
the number of massless fiavors. In the MS scheme, to ~ the thermodynamic potential 
is 
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where G = Gq — 0.536Nf + Nf In Nf, Gq = 10.374 ± .13, and A is the renormalization 
subtraction point. We take the scale dependence of the strong coupling constant, = 



Q;s(A) to three loop order [0. In MS scheme, the thermodynamic potential is manifestly 
gauge invariant. However, there remains a dependence on A/yU that, in principle, can be 
freely chosen. Although the results are significantly sensitive to the particular value of 
this free parameter, this choice is tightly constrained by the physics. (It is not especially 
sensitive to what order in the loop expansion one used for the running of the coupling 
constant, or for the pressure.) 

With this "perturbative" equation of state — used far outside the region of validity for 
perturbation theory — reasonable values turn out to be in the range 2 < A//i < 3. This 
comes about from the most elementary considerations of how to match the quark and 
hadronic equations of state. 
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Figure 1. The total pressure up to order 
~ a^, relative to the pressure of an ideal 

gas, Pfree^ ^ 

function of fi for different 
choices of A//x. 



Figure 2. The total pressure, relative 
to the pressure of an ideal gas, pjree', 
including terms to order ~ and to 
order ~ a^, as a function of /i; A = 2/i. 



In Fig. |I| we show the pressure for some illustrative choices of A//i = 1,2,3. In each 
case, the pressure vanishes at some /x = fic- As the chemical potential of any field must 
be greater than the mass, and as the chemical potential of quarks is one third that of 
baryons, the hadronic pressure vanishes at fimin = 939/3 = 313 MeV (This neglects the 
small effects of nuclear binding, which contribute at most ~ 15/3 = 5 MeV). Then Hc 
must be greater than fimin] from Fig. |1|, this requires that A is less than 3fi. Conversely, 
fj^c cannot be too large. Dense hadrons exert pressure, so before the quark pressure goes 
negative, a transition to massive quarks and hadronic matter occurs. For example, when 
A = n, one can show that hadronic matter must persist to very high densities, maybe 
thirty times nuclear matter density. This is unlikely; indeed, below we see that the worry 
is not how small the hadronic pressure is, but how big. 

Thus, reasonable values of A are in the range of /i — > 3fi. Choosing the value A = 2/i, 
we can get an idea of how well perturbation theory works for large chemical potential. 
At /X = 1 GeV, where as = 0.31, the perturbative corrections to ideality are ~ 30% ; at 
II = 100 GeV, where as = 0.095, they are still 7%. Even so, as we argued in Ref. [|10|, 
the perturbative series at T = and 7^ may be much better behaved than at T 7^ 
and /i = |T^. While for T 7^ 0, yU = we have an expansion in ^/a^, for T = 0, /i 7^ 
it is a perturbative series in as and aslog{as)- Besides, at nonzero density (and T = 0), 
the ~ as, and ~ corrections both act to lower the pressure. It is encouraging that 
the ~ corrections are smaller than those at ~ as down to /i = 0.54 GeV. In contrast, 
at nonzero temperature (and /i = 0), terms oscillate in sign, and the naive perturbative 
series is not well behaved until extremely high temperatures, ~ 10"^ GeV||17||. 

Of course a much better estimate of the quark pressure could be obtained if the terms 
~ al{\og^{as),\og{as),l) were computed. Unlike nonzero temperature — where only 
the ~ a^(log^(Q;s), log(as)) are sensible to compute, for dense QCD all terms in the 
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perturbative expansion of the thermodynamic potential are well defined. 

This does not imply that a given value of a^, which is adequate to compute the thermo- 
dynamic potential, works equally well for all other quantities. In particular, the gaps for 
color superconductivity are nonperturbative, (p ~ exp(— l/^/aJ) and much smaller 



values of appear to be required to reliably compute them pO[. Perturbative results, 
unjustly extrapolated down to ~ 400 MeV, suggest that these gaps are ~ 30 MeV [pi]]. 
As the relative change in the thermodynamic potential is only ~ then, for the 

equation of state in QCD, color superconductivity is a negligible effect. 

The structure of a quark star is determined by the solution to the Tolman-Oppenheimer- 
Volkov (TOV) equations [H- Given the equation of state p = p(e), one can integrate the 
TOV equations from the origin until the pressure vanishes at the edge of the star, obtaining 
a characteristic curve in the mass-radius diagram. In matching the pressure of quarks onto 
that of hadrons, there are two possible scenarios. In the first, the pressure of hadrons 
rise quickly, smoothly matching onto that of quarks. Since the pressure of the hadronic 
and quark phases are similar, so are the densities (as the slopes of the pressure); then the 
chiral phase transition is either weakly first order, or just crossover. In the second case, 
the hadronic pressure is assumed to always remain small. This gives a strongly first order 
chiral transition, as one goes abruptly from a phase of dilute hadrons to one of dense 
quarks ||2^ . A sketch of both cases is depicted in Fig. |^, where the equation of state for 
massless quarks is taken to be the one with A = 2/i. 




Figure 3. Sketch of the total pressure, 
relative to the pressure of an ideal gas, 
Pfree, for a Weak (upper line) and a 
strong (lower line) chiral transition as 
discussed in the text. 



Figure 4. Mass-radius relation for pure 
quark stars (solid lines) and hybrid stars 
(dashed lines) for different critical den- 
sities for the chiral transition. 



Figure | shows the mass-radius diagram for these two cases. A weakly first order chiral 
transition is illustrated by the choice of A = 3fi. The quark pressure vanishes at fic = fJ^min, 
so that the pressure of dense hadrons grows immediately, and becomes equivalent to the 



5 



pressure of dense quarks. The numerical solution of the TOV equations for such an 
equation of state generates hybrid stars with maximum masses and radii compatible with 
those from conventional neutron star models, M ~ 2 Mq and i? ~ 12 km. 

For a strongly first order chiral transition, as with A = 2/i, there is a branch along which 
stars are predominantly composed of dilute hadrons, with small pressure. These look like 
ordinary neutron stars. In addition, there is a new branch in the mass-radius diagram. 
These are predominantly composed of very dense quarks with a thin hadronic mantle. The 
pressure is like that of QCD with a larger bag constant (/imm is larger). Working out how 
the maximum mass and radius scale with an (effective) bag constant, one finds maximum 
masses M ~ 1 Mq and radii i? ~ 6 km. The new branch was suggested previously by 
several groups 

Since we have no plausible hadronic equation of state with small pressure, we simply 
ignore it, and take the quark equation of state down to the point where it vanishes. In 
comparison, in Figure ^ we also show the mass-radius relation for a hadronic equation of 
state (a mean- field equation of state with parameter set TMl, including hyperons PB[]). 
Matching the hadronic equation of state onto the one for the quark phase is done by fixing 
a critical density. For critical densities about twice that of nuclear matter, we find a new 
branch appears, characterized by a second bump in the mass-radius curve. 

Experimentally, Pons et al ^ find that for an isolated, nearby neutron star, the masses 
and radii are exactly what one expects for a new branch of quark stars, with M ^ .QM© 
and R Q km. 



The difficulty is really in understanding how hadrons can have a small pressure 
Akmal, Pandharipande, and Ravenhall have found that to a very good approximation, 
in pure neutron matter the energy per baryon (minus its rest mass) is approximately linear 
in the baryon density, n: 

m = m=15 MeV ( — ) . (3) 



A n V^O' 

where e is the energy density and Uq ~ .16 baryons/fm^ is the saturation density for 
nuclear matter. These authors used a wide variety of different models for the nucleon- 
nucleon interaction, and find that up to n ^ no, all models are in close agreement. Even 
at such "low" densities, this behavior itself is already signs of a highly nonideal Fermi 
liquid, since free fermions give e/n — m ~ n^/^, not ~ n. By two times nuclear matter 
density, the discrepancies between different models is already ~ 50%, and grow sharply 
with density. 

The energy on the right hand side of (|^), ~ 15 MeV, is very small on hadronic scales. 
It is much smaller than the pion decay constant, which one might expect is the natural 
scale. What is remarkable about the results of Akmal et al. E^ is that the energies are 



small not only for nuclear matter, which has a nonzero binding energy, but even for pure 
neutron matter, which is unbound. This suggests to us that the smallness of the right 
hand side of eq. (^ is a generic property of baryons interacting with pions, etc., and is 
not due to any special tuning. Using this result, however, we can reliably compute the 
pressure of pure neutron matter at nuclear matter density: 

d / € \ Ti^ / Tl \ 

Phadron = TT " = McV — = .04 pjree — (4) 

on \nj no \noJ 
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Here p free is the pressure for three flavors of massless quarks; for the purposes of discussion, 
in pfree we have taken = 313 MeV. Thus at nuclear matter density, the pressure is 4% 
that for ideal quarks, which is comfortably small. For the hadronic pressure to remain 
small relative to the quark pressure, the transition must occur at remarkably low densities, 
maybe twice that of nuclear matter. 

We do not have a good reason for why the hadronic pressure should remain small to 
the chiral transition. But since we don't understand why it is small in nuclear matter, we 
feel free to extrapolate; maybe the hadronic pressure is suppressed in the limit of a large 
number of colors. Such "implausible" equations of state, with a small hadronic pressure, 
can be constructed mathematically ^6 . 



If small and dense quark stars occur, then their production is probably different from 
those of ordinary neutron stars. Models of supernova collapses suggest that the mass of a 
proto-neutron star formed by the collapse of the iron core is in excess of IMq as the initial 
mass of the iron core is well above IMq. Neutron stars can not collapse to quark stars, as 
they have to expel most of their mass to achieve a stable quark star. White dwarfs do not 
collapse to more compact objects in type la supernova as the explosion destroys the whole 
star. In the early universe, blobs of strange quark matter formed during the quark-hadron 
phase transition can not be larger than 10~^Mq due to the horizon at t = 10~^s after the 
big-bang . 



Nevertheless, there might be some exotic ways to produce quark stars. If white dwarfs 
can develop an iron core by some (so far not established) reason, the collapse of this core 
can form a low-mass quark star. Alternately, if inflation occurs at the quark-hadron phase 
transition, then larger strange quark stars, with M = 0.01 — IOMq, can be formed p8| . 



The existence of a new class of small and dense compact stars is an exciting possibility. 
While superconducting condensates do not matter for the pressure, they do dominate 
their electromagnetic response. Whatever the magnitude of the gaps for color supercon- 
ductivity, versus those for hadrons, there should be a multitude of ways in which stars 
composed of color-flavor locked quarks differ from those of superconducting (and super- 
fluid) baryons. After all, in this scenario the stars are never really hybrids, but either 
mainly hadronic, or mainly quark. 
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ment of Energy for their support under Contract No. DE-AC02-98CII10886; E.S.F., for 
the support of CNPq (Brazil); and J.S.B., for the support of RIKEN and BNL. 
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